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Abstract 

Aspects of the phase change of the two-level pairing model are investigated in 
the semi-classical treatment by using the variational approch with the mixed-mode 
coherent state. In the classical limit, h — >• 0, the sharp phase transition appears and 
the two phases exist in the region where the force strength is larger than a certain 
critical value. However, it is shown that, in the semi-classical treatment, the above- 
mentioned behavior of the phase change disappears in the analytical and the numerical 
treatments. As a result, new understanding about aspects of the phase change in this 
model is given. 
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§1. Introduction 



The phase structure and the phase transitions between various phases in condensed mat- 
ter systems, such as nuclear and/or quark matter as well as the condensed matter in atoms 
and molecules, are important concepts to understand the nature of the matter under con- 
sideration. In the infinite systems, the clear phase separation is often realized and, in that 
case, the sharp phase transition may occur among the different phases. On the other hand, 
in finite systems such as nucleus, the phase change is not so clear. This situation can be 
demonstrated 1 ^ in the exactly-solvable many-particle models such as the Lipkin model 2 ) or 
the pairing model. 3 ) However, the mean field approximation in the time-independent the- 
ory, which often leads to the classical treatment, gives the fictitious sharp phase transition. 
Addition to the above mentioned problem, the two phases, namely, the condensed and the 
non-condensed phases, appear in both the Lipkin and the pairing models in the classical 
treatment, if the force strength is larger than a certain critical value at zero temperature. 

In this paper, an exactly-solvable two-level pairing model is treated to understand the 
nature of the phase change and the existence of two phases for the force strength at zero 
temperature in the classical limit. The two-level pairing model is regarded as the su(2) ® 
sw(2)-algebraic model. For the original many-fermion system, the trial state for the variation 
is restricted to the BCS state which gives the mean field approximation and the unique 
classical counterpart for the original quantal system. In order to avoid the usual mean 
field approximation of the BCS theory, the two-level pairing model can be expressed in 
the four kinds of boson operators by using the Schwinger boson representation of the su(2)- 
algebra. In this boson representation, the various trial states are possible to obtain a classical 
counterpart. 4 ) Thus, in this paper, a certain wave packet which we call the mixed-mode 
coherent state 5 ) is adopted to the variational calculation. This state gives the semi-classical 
approximation which retains H 

This paper is organized as follows: In the next section, the two-level pairing model is 
introduced by means of the Schwinger boson representation. In §3, the classical counterpart 
for the two-level pairing model is given by using the mixed-mode coherent state and the 
basic equations of motion in this approach are also derived. In §4, the energy minimum 
point is investigated in the classical limit h — > and in the semi-classical approach with 
finite value of h. Then, it is shown that the sharp phase transition does not occur with 
finite h, while the sharp phase transition occurs and two phases exist in the certain region in 
the classical limit. The method to obtain an approximate solution for the energy minimum 
is given analytically in §5 together with appendix A. The behavior of the phase change is 
explained with unknown concept in §6 and the discussion about the ground state energy is 
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given in §7. The last section is devoted to a summary. 

§2. Two-level pairing model in many-fermion system and its counterpart in 

the Schwinger boson representation 

The model discussed in this paper is two-level pairing model in many-fermion system. 
The two levels are specified by a = + (the upper) and o = — (the lower), respectively. 
The difference of the single-particle energies between the two levels and the strength of the 
pairing interaction are denoted as He (> 0) and h 2 G (> 0), respectively. The Hamiltonian, 
H, is expressed in the form 

H = e(S (+) - 5bH) - G(S+(+) + S+(-))(S-{+) + S-H) • (2-1) 

Here, S (a), S + (a) and S-(cr) are defined as 

S (a) = (H/2)(N(a) - Q a ) , (2-2a) 
S + (a) = (ft/2)P» , S» = (h/2)P{a) , (2-2b) 

N(a) = > n(a) = (2j a + l)/2 , (2-3a) 

m=-j a 

f>) = E ^(^)(-) J ' CT " m 5- m (^) , H°) = E (-) J ^ m c_ m ( ( x)g m ( ( T) . (2-3b) 

m=—j a m=—j <7 

The operator (c^(cr), c m (cr)) denotes fermion operator in the single-particle state (a, j a ,m;rn - 
—jo-, —j a + l, ■ ■ ■ , j a — 1, j a ). Of course, jo- is a half-integer. Clearly, iV(cr) and (P*(a),P(a)) 
denote the fermion number and the fermion pair operators in the state a, respectively. 

We know that the set (S± ; o(a)) obeys the sw(2)-algebra, and then, essentially, the pre- 
sented model is governed by the su{2) (g> sw(2)-algebra. Clearly, the addition of the two sets 
also obeys the sw(2)-algebra: 

3 ± ,o = 5 ±> o(+)+£t >0 (-) • (24) 

The Hamiltonian (|2T|) has three constants of motion. This can be shown from the following 
relation: 

[ S\+) ,H] = [ S\-) , H ] = [ S , H ] = . (2-5) 

~2 ~ 

Here, S (a) denotes the Casimir operator for (S± ) o(cr)). It should be noted that the Casimir 
operators for (S± j o), S , does not commute with H: 

[ S 2 , H } ? . (2-6) 
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A general framework of the su(2) <8> sw(2)-algebra can be formulated in a form of the 
Schwinger boson representation, in which four kinds of boson operators are used. A possible 
example of this formulation was performed by the present authors (J. P. and M.Y.) with 
Kuriyama. 5 ** Hereafter, we will refer it as (A). In (A), the su(2) ® sw(2)-algebra is formu- 
lated in terms of the bosons (a* , a a , 6* , b a ; a = ±). In this framework, we can construct a 
counterpart of the model presented in this section in the Schwinger boson representation. 
The sw(2)-generators (S± j o(a); a = ±) correspond to S± t o(cr) in the form 

S± >0 (a) — > S±, (<r) , (2-7) 
S,{a) = {h/2){ala u -blb u ) , 

S + (a) = ha% , S-{a) = hb* a a a . (2-8) 
It should be noted that the notations are different from those in (A). Of course, we have 

S = (H/2)(a* + a+ - b\b + + a*_a- - b*_L) , 
S + = h{a* + b + + a*_b-) , 

S_ = h{b\a + + b*_a-) . (2-9) 

Associating with the above operators, we can define the following operators: 

S(a) = (H/2)(ala a + b*X) . (2-10) 

This operator satisfies, for the Casimir operator S (a), the relation 

S\a) = S(a)(S(a) + h) . (2-11) 

In the case where the seniorities are equal to in the fermion system, we should set up the 
correspondence 

(h/2)f2(a) — ► S{a) . (2-12) 

In the Schwinger boson representation, there does not exist the concept of the fermion 
number explicitly. The definition of 5*0(0") given in the relation (j2-2aj) and the correspondence 
(|2-12j) permit us to introduce the fermion number in the Schwinger boson representation in 
the form 

HN(a) = 20(a) + S Q (a)) = 2Ha* a a a . (2-13) 

Therefore, the total number is given as 

HN = h(N(+) + A>(-)) = 2(S(+) + S(-) + S ) 

= 2h{a\a + + a*_a-) . (2-14) 

From the above form, we can see that in the present case for both levels, even fermion 
numbers should be occupied. In the case of odd fermion numbers, we must take into account 
the case where the seniorities do not vanish. 
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§3. Classical counterpart of the su(2) (g) sw(2)-algebra 



We continue the recapitulation of (A). In (A), we showed that the orthogonal set for the 
su(2) <S> sw(2)-algebra is obtained with the use of the set (R±,o,T± ;0 ) as a supporting role. 
The set is defined in the form 

R = (h/2)(a* + a + -a*a_) , 

R + = h(a* + a_ + S;S_) , R- = h(a*_a + + b*_b+) , (3-1) 
f = (h/2)(a\a + + S* 5_ + a*La- + b\b + + 2) , 

f + = h{a* + b*_ - a*_b* + ) , f_ = h(b_a + - b + a_) , (3-2) 

The sets (R±,o) and (2±,o) obey the su(2)- and the su(l, l)-algebra, respectively, and it is 
characteristic that three sets (5±,o), (R±,o) and (T±,o) commute with one another. Further, 
the Casimir operators for the three forms are identically equal to one another. 

In (A), we defined a wave packet, which we called a mixed-mode coherent state, as follows: 

|c) = t/r'exp (-\W\ 2 /h) ■ exp ({V t /hU t )f + ^ • exp {{V s /hU s )S^j 

x exp (jy r /nU r )R+ > j ■ exp {[\fh/2WU r U s /U^j b*_^j |0) . (3-3) 

The state |c) is normalized and U r , U s and U t are real and W , V r , V s and V t complex. They 
obey the condition 

U? - \V t \ 2 = 1 , U 2 S + \V S \ 2 = 1 , Ul + \V r \ 2 = 1 . (3-4) 
We can prove the following relation: 

a + = U t A + + V t B*_ + ^(WV r V s U t + WE/^) , 



a_ = U t A_ - V t B\ + \/-(WU r V s U t - W*V r *U s V t ) , 



b + = U t B + - V t A*_ + ] Jl(WV r U s Ut - w*u r v;v t ) , 

b_ = U t B_ + V t A\ + ^(WU r U s Ut + W*V r *V s *V t ) . (3-5) 
Here, (A a , A*) and (B a , B*) denote boson operators satisfying 

A a \c) = , B a \c) = . (3-6) 
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As a possible parameterization, in (A), the following form were presented: 
m _ v /^ e -j(V'-^o+V'(-))/2 

K = y/(S + S(+) - S{-))/2Se- iW+) -^- ))/2 , 



V s = ^{S + S Q )/2Se-^ , 

V t = ^/(S(+) + S(-) - S)/2(S + h)e~ i{ ^ +)+ ^- ))/2 



U r = y/(S-S(+) + S(-))/2S , U s = V(S-S )/2S 



U t = y/(S(+) + S(-) + S + 2h)/2{S + h) . (3-7) 

In (A), we prove that (ip(cr), S(a)), (tpo, S ) and (tp, S) are in the relations of the angle and 
the action variables in classical mechanics. Further, we have the following relation: 

S>0, S<S <S, \S{+)-3{-)\ <S< 3(+) + S(-) . (3-8) 

The above is nothing but the relation that implies the coupling of two sw(2)-spins. 

With the aid of the relation fj3-5j) . we can calculate the expectation value of any operator 
composed of the present bosons with respect to |c): We denote (c|0|c) = (0) c . For example, 
we have 

(S(a)) c = S(a) , (3-9) 
(So)c = So , 



(S±) c = ^S 2 - S 2 e^ . (3-10) 

The relation (|3-1(J|) is identical with the classical counterpart of the Holstein-Primakoff boson 
representation, and then, including the relation ( |3-8jl . our present formalism may be called 
the classical counterpart of the su(2)(g>sM(2)-algebra. Concerning S+S-, the most interesting 
point is as following: 

(S + S^) c = ((£+(+) + 5+(-))(5_(+) + 5_(-))) c = S 2 - S 2 + h(S + So) . (3-11) 
(S + ) C (S_) C = S 2 - S 2 . (3-12) 

We can see that in the relation (j3-ll|) . the quantal fluctuation is exactly taken into account. 
Further, we have 

(S (+) - S (-)) c 
(S(+)-S(-))(S(+)+S(-) + h)S Q 

s(s + h) 

^(S - S(+) + S(-) ) (S+ S(+) - S(-) ) (S(+) + S(-) - S) (S(+) + S(-) + S+ 2h) 
+ W^h) W ' 

(3-13) 
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Under the above preparation, we are able to obtain the classical Hamiltonian of H which 
is the counterpart of H shown in the relation (|2-1|) : 

H = e{S (+) - S (-)) - G{S + (+) + £+(-))(£_(+) + &_(-)) . (3-14) 

The Hamiltonian H obeys the condition 

[ S(+) , H } = [ S(-) , H] = [§ , H] = . (3-15) 

The relation (|3T5jl corresponds to the relation (|2-5jl . It is enough to calculate the expectation 
value of H for |c): (IT) C = -f/". In this case, the forms (j3-ll|) and ()3-13j) are useful. Of course, 
the classical expression of the fermion number is given in the form 

hN = h(N) c = 2(S(+) + S(-) + S ) . (3-16) 

For the above, the relation (j2-14j) is used. Corresponding to the Hamiltonian H, H satisfies 

[ S{+) , H] P =[ S(-) , H } P = [ So , H }p = . (3-17) 

Here, [ , }p denotes the Poisson bracket. As was already mentioned, the parameters 
(ip(a), S(a)), (x/jq, So) and (if), S) are in the relation of the canonical variables in classical 
mechanics. From the relation (|3-17|) . we can see that <S'(+), S(-) and So are constants of 
motion and the time-evolution of the angles and ipo are given in terms of 

i>(+) = [1>(+), H] P , ^(-) = [^(-) , H] P , iPo = l*Po,H}p. (3-18) 

The variables (ip, S) satisfy the Hamilton equations of motion: 

j> = [if> , H] P , S=[S,H] P . (3-19) 

In order to avoid unnecessary complication, we treat the following case: 

s(+) = s(-) (=l) , L = \n, (o(+) = o(-) = n) 

So = . (3-20) 

In this case, the relation (j3-16J) is reduced to 

hN = AL , i.e., N = 2(2. (3-21) 

The relation (|3-21j) shows us that we are interested in the closed shell system. Under the 
expectation values (|3T2j) and ()3T3|) and the condition (|3-20p . H can be reduced to 



H = e[l 



( 1 - — ^— ] J(2L -S)(2L + S + 2h) cos^ - G(S 2 + hS) . (3-22) 
\ S + nj 



7 



The Hamiltonian ()3-22j) can be rewritten as 

H = e(l-^j V(2A) 2 - o- 2 cos V - Ga{cr - h) . (3-23) 
Here, A and a are defined by 

A = L + h/2 , a = S + ft. (3-24) 

Concerning the parameters A and a, the following relation should be noted: Since there exists 
the restriction \S(+) —S(—)\ < S < S(+) +S(—) classically and quantum mechanically, for 
the case S(+) = S(-), we have 

h<a<2\. (3-25) 
Since (ip, a) is also in the canonical relation, we have the following Hamilton's equations: 



i- = [ i, , H } P = = € 



cosip 



v 2 V a) v/(2A) 2 - a 2 

-G(2a - K) , (3-26a) 

o = \ a , H] P = -^- = e( 1-- ] v /(2A) 2 - ( r 2 sin^ . (3-26b) 



§4. Determination of the energy minimum point 

One of our interests in this paper is to investigate how the energy minimum point induced 
by the mixed-mode coherent state |c) changes in the phase space (if), a) as a function (e, G). 
The energy minimum point should be stationary in the phase space, and then, we have the 
condition 

^ = 0, (7 = 0. (4-1) 

The relation (j3-26b|) gives us 

smip = , i.e., cost/? = — 1 or + 1. (4-2) 

It may be self-evident that the form (j3-23J) supports 

cosV> = -l. (4-3) 

The angle if) does not depend on (e, G). Under the conditions (|4-lj) and ()4-3j) . the relation 
()3-26a|) is reduced, after some calculations, to 

a 3 (e 2 -4G 2 ((2\) 2 -a 2 )) 
=H (e(2A) 2 - GaV(2A) 2 - d 2 ) (e + 2G V / (2A) 2 - a 2 ) . (4-4) 
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By solving Eq. (j4-4j) . we can determine the value of a which makes the energy minimum. As 
a possible solution of the latter of Eq. (|4-l|) . we obtain a = 2A, which is independent of ip 
and e. Substituting a = 2 A into the former of Eq. (j4-l|) . we have cos i/j = — (4GA/e)(4A — 
K)/{A\ — 2h) ■ \/\5. Here, 5 denotes an infinitesimal parameter defined as a = 2A — 5. Then, 
under 5 = 0, cosip = 0, i.e., ip = ±vr/2, which is independent of e and G. But, this solution 
cannot connect to the case o < 2A. Therefore, we do not adopt it. If we adopt the solution 
cos ip = — 1, clearly, 4GX/e becomes infinite under the condition 5^0. This means that 
our two levels are degenerate. 

Let us search the solution of Eq. (j4-4|) . This relation can be simplified in the form 



p 3 (i - x 2 (i - p 2 )) =* (i - x ~y^?) (i + ^ 



P 2 ) ■ (4-5) 



Here, fc, p and x are defined as 

k = h/2X , (4-6) 
p = <x/2A , x = AGX/e . (4-7) 

The relation (jH-24j) and the condition L > give us 

< k < 1 . (4-8) 

Especially, we have 

fe — > , if A — > oo , i.e., L — > oo , 

h 

k = 1 , if A = - , i.e., L = . (4-9) 

If we take into account the correspondence to the original fermion system, the minimum 
value of i? is equal to 1, and then, L = h/2 and we have 

< k < 1/2 . (4-10) 

However, formally, we have the minimum value of L = 0. Further, the relation ()3-25J) gives 
us 

k<p<l. (4-n) 

Of course, we have 

x > . (4-12) 

First of all, on the basis of the relation (J4-5)) . we treat the case h (the Planck constant) 
= 0. In this case, we have a simple solution: 

any value (> 0) , for p = , 
x = { (4-13) 
, , for < p < 1 . (x — > oo , if p — > 1) 
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Fig. 1. The variable x is depicted as a function of p with h = 0. 



Inversely, the solution (|4-13jl are given as 



P 



, for any value of x(> 0) 
■ , for x > 1 . 



v" 



(4-14) 



In Figd we see a straight line consisting of solid and dotted lines along x-axis. This is one 
phase which is conventionally called as normal conducting phase. Further, we see that at the 
point i= Ion the straight line (x-axis), a new phase is born toward the rectangular direction 
and it grows up to a branch, which is conventionally called as superconducting phase. In 
FiglU this phase is specified by a dot-dashed curve. The above is a conventionally accepted 
fact. In the region x > 1, there exist two phases and the energy of the superconducting 
phase is smaller than that of the normal conducting phase. From this reason, we chase the 
system under investigation along the solid to the broken line, the connection of which is not 
smooth. In this sense, conventionally we call this situation the sharp phase transition. 

Our next task is to investigate the system under the condition h =finite. For this purpose, 
it may be enough to solve the quadratic equation with respect to x, (|4-5|) . In order to see 
the whole feature of the solution, for a moment, we treat the variables p and x in the region 



-1 < P < 1 , -oo < x < oo . (4-15) 
Two solutions of Eq. (j4-5|) for the case fc ^ and 1 are given exactly in simple forms: 



x 



l p 3 - k 



1-p 2 p* [P 2 



(4-16a) 
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Fig. 2. The variable x is depicted as a function of p with h ^ (fe = 0.75). 



1 



.r 



(4-16b) 



The value of p at the cross point of the solutions ()4T6a|) and (J4T6bj) . which we denote pp, 
is given by 

1 / 7„ \ 



PP 



1_ 

51 



1 + 



51 



+ 



54 (2) 



V 



V 



1 + JL h 
L T 54 l 2 



(4-17) 



Figure El shows the behavior of x as a function of p. The solid curves show the behaviors 
of x given in the relation (j4-l(iaj) in the region (0 < p < 1 , 0<a;< 00) and in the other 
region, respectively. The dot-dashed curve shows the behavior of the solution (|4T6b)l . We 
can see that between Figs|T] and El there exists quite big difference. In §6, we will discuss 
this difference. At first sight, FigsnandOltell us that, under the limit h (fc) — > 0, the solid 
curve in the region fc 1 / 3 < p < 1 of FigJ21 approaches to the solid line and the broken curve in 
FigC] and the solid curve in the region < p < fc 1 ^ 3 is compressed and becomes the dotted 
line in FigO Through the discussion in §6, the above view will be answered in the negative. 
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§5. Approximate solution for the relation determining the energy minimum 

point 

As was shown in Eqs. (j4-16a|) and (j4-16b|) . we can express x in terms of p. However, in 
order to understand the behavior of the system, for example, the change of the energy, in 
relation to the interaction strength, we must express p as a function of x. Numerically, it is 
possible and easy But, analytical expression may be also helpful for the understanding of 
the system. In the present case, the exact form cannot be almost expected and we will try 
to the approximate form. 

Let us investigate the case of the smooth curve in the region fc/2 < p < 1 in Fig|21 First, 
we note Eq. (|4-16a|) . This equation can be rewritten in the form 

P - fe 1/3 fl^p- p 2 + fe 1/3 p + fe 2/3 f . 

X (l-p)(p-fc/2)'Vl + p' P 2 ' 1 j 

Then, we define the quantity y in the form 



V= ^- JTW^ ' X - (5 ' 2) 

With the quantity y, the relation (|5Tjl is formally reduced to the quadratic equation for p: 

V = f{p) ' f{p) = (1 -p)(p- */2) ' (5 ' 3) 
Formally, from the relation ()5-3|) . we have 

P = F(y) , 

F(y) = % " 2 2feV3 • (5-4) 

(i + 1) y - 1 - y (i - 1) y 2 - 2 (i + 1 - 2fc x / 3 ) y + 1 

If x — » — oo, and +oo, i.e., y — > — oo, and +oo, we have the relations F(— oo) = fe/2, 
fc 1//3 and 1. Therefore, at the three characteristic points, we are able to obtain the exact 
results. From this reason, we do not adopt another branch for the solutions of the quadratic 
equation. 

Next, we consider the case of the points except x — > — oo, and +oo. First, we define a 
series {p^°\ p^\ p^ 2 \ ■ ■ ■} obeying the following recursion relation: 

k = 1,2,3,- ■■ . (5-5) 
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Here, we note that the function f(y) is monotone-increasing. Then, if there exist the limiting 
value limfc^oo p^ = p, p is determined as a function of p by the condition 

p = F(y) , y = f(p) . (5-6) 

The relation ()5-6|) is nothing but the set of the relations (|5-3jl and ()5-4|) . Under appropriate 
choice for the initial value p(°\ iteratively, we obtain p^ from k — I and after infinite 
iteration, we may obtain the exact solution which do not depend on p^ . If we stop the 
iteration half way, the result is an approximated one which depends on p(°\ As a possible 
choice of the initial condition, we adopt the value of p satisfying the form (j4-16a|) at x = 1, 
which will be denoted as p c : 

1 Pl ~ k (5-7) 



^rr^ c pl( Pc -k/2) 

The relation ()5-7p leads us to 



P 5 



Since the right-hand side of the relation ()5-8|) is increasing, and then, k and p c are in the 
relation of one-to-one correspondence. Therefore, instead of k, we can use p c defined in the 
relation (|5-8|) . Therefore, the cases x = — oo, 0, 1 and +oo do not depend on the iteration. 

§6. Examination in the case fe — > and 1 and its related phase transition 

In §4, we showed the behavior of x with respect to p in the case where k does not approach 
to and 1. In this section, we discuss the case where k approaches to and 1. First, we 
treat the case k — - > 0. In the case p < p c , we have 

x= H —, ^^^^=, (fe->0). (6-1) 

Further, in the case p = p c , we have 

v/wfp^ (p c - 1) vw! 

namely, 

p c — > . (fc -> 0) (6-2) 
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Fig. 3. The variable x is depicted as a function of p with — » (k = 0.3, 0.2, 0.1). 



In the case where p is a little bit larger than fc 1 / 3 , we can put p = (ak) 1 / 3 (a > 1). Then, 
we have 

ak — k 



x 



y/l - (ak) 2 ' 3 (ak) 2 / 3 ((ak) 1 / 3 - 
a — 1 

v/l - (afe) 2 / 3 (a - |a 2 / 3 fe 2/3 ) 
The above relation leads us to 

1 



1--, (fc-0). 



p^0 



;r 



1 - 



< 1 - - < 1 



(fe^O) . 



(6-3a) 



(6-3b) 



The relation ()6-3b|) tells that under the limit k — > 0, p — > and a; takes arbitrary value 
between and 1. From the above three cases, we see that x in the region fc 1 / 3 < p < 1 in 
Figgis reduced to the solid line and the dot-dashed curve in Fig|T] Next, we treat the case 
where p is a little bit smaller than k 1 ^. In this case, we can set up p = (k/a) 1 / 3 (a > 1). 
Then, in the same process, we have 



x 



1 — a 



(k -> 0) . (6-4a) 

The above means the following situation: 

p — ► , x — ► 1 - a , (1 - a < 0) . (k -> 0) (64b) 

The above situation can be interpreted as that under k — > 0, p —> and x takes negative 
arbitrary value and in Figd the straight line between and 1 continues to the negative region 
of x. In the case where p is in the region < p < k/2, we can put p = a(fc/2), (0 < a < 1). 
Then, we have 

(ak/2) 3 - k 



x 



y/l - (ak/2) 2 (ak/2) 2 (ak/2 - k/2) 
8/k 2 - a 3 



\J\ - a 2 k 2 /A a(l - a) 



+oo . (£-►()) 



(6-5) 
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Fig. 4. The variable x is depicted as a function of p with fc — * 1 (fc = 0.8, 0.9, 0.99). 



This situation is very interesting. The solid curve in the region < p < k/2 in FigJ2] 
disappears to +oo. Therefore, the dotted line in FigH] does not correspond to this curve. 
Next, in the case where p is a little bit smaller than 0, we can set up p = — (ak) 1 ^ 3 (a > 0). 
Then, x can be expressed in the form 

ok + k 

x 



y/l - (akf/ 3 a 2 l 3 k 2/3 (a l / 3 k l/3 + |j 



n + j ■*! + -, (fe->0) . (6-6a) 



y/l - (ak) 2 / 3 a + \a 2 l 3 k 2/3 a 
The above means 

p — >0, x — >1 + -, (1 + ->1J . (fe->0) (6-6b) 

a \ a J 

In this case, under the limit k — > 0, p — > and x is positive arbitrary and this situation 
corresponds to the dotted line starting at x — 1 in Fig^ However, this situation comes from 
the negative region of p, which is physically unaccepted. Further, in the case p <C — (ak) 1 / 3 , 



X^l/y/l-p 2 (k^0). 

Next, we treat the case k — > 1. In this case, we divide two situations. One is the case 
k 1 ^ 3 < p. If we put k = 1 — Sk and p = 1 — 8p, we have Sp < 5k/3. Then, we put Sp = aSk/3 
(0 < a < 1). Under the condition that 5k (Sp) is sufficiently small, x can be expressed as 

x = ( 2~ x ^ Sk i — » 3 (r ~ 1 ) > fr" 1>0 V (*-°) ( 6 ' 7 ) 



In this case, under the limit — > 1, p — > 1 and x is positive arbitrary. The above means that 
the solid curve between p = k 1 ^ 3 and 1 becomes the straight line on p = 1. On the other 
hand, for the region k/2 < p < k 1 ^ 3 , the limit k — > 1 gives 



x — t/irf^W- (6 ' 8) 
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The above shows that the solid curve between p = k/2 and fc 1//3 in Fig|2] becomes the curve 
showing by the relation (|6-8|) . The other cases do not have characteristic change. However, 
as was already mentioned, k — *■ 1 corresponds to L — > 0. In the system in which the seniority 
is zero, L starts from L = 1/2, and then, the present case is not so interesting. But, in the 
case where the seniority is not zero, it may be interesting. 

On the basis of the above-shown results, let us consider the phase transition in our model. 
In the present case, physically interesting feature can be seen in the solid curve between 
p = fc 1//3 and 1 in FigJ21 We cannot find any non-smoothness or discontinuity. Of course, 
under the limit k — > 0, the result appears as the phase transition which is conventionally 
accepted. Therefore, in the finite h and A, we should not expect the phase transition. 
Further, we cannot find any curve which leads us to the dotted line in FigHfrom the side 
of the physically acceptable region. Therefore, the dotted line in FigQ should be negative 
understood. From the above examination, we can learn that the limiting process H — > 
must be carefully performed. All do not follow this process. 

Let our model accepted in the region x < (the repulsive pairing interaction). In FigEl 
we see that there exists a cross point of two different curves. Noting the relation (|4T1J) 
(k < p < 1), we examine characteristic feature of the crossing. There exist two types of 
the crossing which show in Figl^a) andEfb). Both figures are characterized by pp > k and 
k/2 < pp < k, respectively. Here, pp is shown in the relation (|4-17|) . From Fig|SJ we learn 
that there exist three types of non-smoothness and discontinuity, which are illustrated by 

Fig. 6. The three types of non-smoothness and discontinuity for the crossing point pp are illus- 
trated in the case k < pp. 



Fig. 7. A type of non-smoothness and discontinuity for the crossing point pp is illustrated in the 
case k > pp. 
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Fig® Following the change of the interaction strength, the system moves along the solid 
curve. At the present, we do not judge which should be chose in the three cases. In the case 
of Figl^b), we can draw the picture shown in Fig|7| Discontinuity appears in the interaction 
strength. Thus, we can conclude that in our present model, the phase transition appears 
in the case of the negative interaction strength. However, at the present stage, we cannot 
judge if the above phase transition has its realistic meaning. 

§7. Discussion 

Finally, in this section, we discuss some characteristic points of the energy E. The energy 
E is obtained by substituting the solution of Eq. (j4-lJ) into H shown in the relation ()3-22j) : 



E = -e • 2A 



(7-1) 



In FigJHl we show numerical results. In the case x —>■ oo (two levels are degenerate), we have 
p = 1 and E is reduced to 



Ei = -G ■ 2L(2L + h) = -h 2 GQ{Q + 1) 



(7-2) 



The form (|7-2j) is the same as that shown exactly. On the other hand, in the case x = (free 
fermions), we have p = fe 1//3 and E is written as 



E 



e • 2A 1 - fe* 



3 

2 \ 2 




If Q is large, the leading term of E is given by 

he 

E « 2J? 

2 



h 



(J? + 1)3 - 1 



2L + h 



- h i N 




(7-3) 



(74) 



The above form tells us that the lower level is fully occupied. Therefore, we can conclude 
that in the case of large f2, our result is reduced to the exact one and for the problem of 
k — > 0, it may be favorable. 



Fig. 8. The ground state energy normalized by e is depicted as a function of the interaction strength 
x = 4GA/e with k = 1/3 (solid curve) together with Eq and E\. 
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§8. Summary 



In this paper, we analyzed the aspect of phase change in the two-level pairing model 
governed by the su(2) <g> sit(2)-algebra. This model can be described in terms of four kinds of 
boson operators by means of the Schwinger boson representation of the sw(2)-algebras. We 
used the variational approach with the mixed-mode coherent state introduced in Ref.5). This 
state is constructed, as is similar to the usual coherent state, by using the raising operators 
of the su(2)- and the su(l, l)-algebras whose generators are defined by the Schwinger bosons 
used in this model. 

In the classical limit H — > 0, in which the Planck constant is taken to be at the first 
time, if the force strength is smaller than a certain critical value, the trivial phase, which 
is called the non-condensed phase, is realized. At the critical point with respect to the 
force strength, a new branch appears and this branch reveals the existence of the so-called 
condensed phase. However, the non-condensed phase still exists as a result of the trivial 
solution if the force strength is larger than the critical value. This situation is well known 
in the usual Hartree-Bogoliubov approximation and the sharp phase transition is realized in 
this treatment. 

However, this behavior is not true because the sharp phase transition is not realized in 
this many-body system and this fact can easily be verified in this exactly solvable algebraic 
model. Figures El shows the aspect of phase structure. Furthermore, in FigJ31 the aspect of 
the phase change in the semi-classical limit is shown clearly. It is seen in Fig |2 that the two 
branches exist in the physical region in < p < 1, where is included when h — > 0. If h 
approaches to 0, one branch, which is located in < p < fc/2, goes to infinity of x and this 
branch disappears in the limit h —>■ 0. Another branch, which is located in fc/2 < p < 1 
and x > 0, is reduced to the physical branch in FigJTJ which is represented by solid and 
dash-dotted curves. However, the branch of "the non-condensed phase" in the classical 
limit, which is represented by the dotted line in Fig.l in x > 1, appears from the unphysical 
region in p < in Fig|21 So, this solution is not a physical solution originally. Thus, this 
"non-condensed phase" is not realized in the physical system. Therefore, the sharp phase 
transition does not occur if h is taken as finite value. 

In our model, there exist a cross point of two branches represented by two curves in Fig|21 
in the negative region of x. Because the variable x is proportional to the force strength G, 
this cross point is realized in the case of the "repulsive pairing". Thus, the sharp phase 
transition may be realized in the case of the negative interaction strength. 
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Appendix A 

A method for solving the relation A5-5}) 



Let an infinite series {y(°\ y^\ y^ 2 \ • • • } obey the following condition characterized by a 
function f(x, y): 

y« = (f(x,y^)Y . (* = 1,2,3,.-.) (A-l) 

Here, v denotes real parameter and x plays a role of real variable (0 < x < oo). Further, we 
regard f(x,y) as a monotone-increasing or -decreasing function. If there exists the limiting 
value lim^oo ?/ & ) (= y), y is determined as a function of x by the condition 

y = (f(x,y)T . (A-2) 

The above means that the variable y governed by the relation (|A-2|) is expressed as a function 
of x by calculating the limiting value lim^oo t/" of the series 
{y(° \ y^\ y( 2 \ ■ • • } obeying the condition (|AT[) . On the basis of the above idea, in this 
Appendix, we develop a systematic method for obtaining an approximate expression of y 
as a function of x. Of course, x and y obey the relation (jA-2jl . and hereafter, we treat the 
following case: Through the relation (jA-2|) . we are able to obtain the values of y at x = 
and at x — > oo. These two values are denoted as y and y^, respectively: 

2/o = (f(0,yo)T , 2/oc = (/(oo,?/ 00 )) ly . (A-3) 

Hereafter, we denote any quantity X at x = and x — > oo as X and Aoo, respectively. The 
derivative (dy/dx)o is given as 

dy\ _ (|£)o ( A4 ) 



»y \ 9 yJo 

The quantities (df/dx)o and (df/dy)o are obtained by putting x = and y = yo for df/dx 
and df /dy as functions of x and y, respectively. 
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Let us set up the following form for y^: 

w = (v&>* + y5<*Y (A . 5) 

Independently from a, we have 

y { o ] = Vo , y£> = Voo ■ (A-6) 
We choose a so as to make y^ satisfy the relation (jA-4|) : 



a = uy 



Voo 

Vo 



1 



(^)o 

'df 



\dxJo 



(A-7) 



Therefore, naturally, we have 



\ dx ) \ rfx 



(A-8) 

\"""/ Q 

Through the relation (|A-lj) . we obtain y^ 1 ' from t/ ) given by the form (|A-5j) . By performing 
the above process successively, we can calculate y^ k > for any k. If k is larger, the accuracy of 
the approximation may be better. We can prove the following form: 



Jk) _ (0) _ (*) _ (0) _ 

(fl) = (^) = (£) . (, = 2 , M ,...) 
\ ax / Q \ ax J \ax J 



(A-9) 



Therefore, we can conclude that, for any fc, the behaviors of i/ fe ) near x = and x — > oo do 
not change from those for fc = 0. 
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